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A theory of magnetic short-range order for high-T c cuprates is presented on the basis of the one-band 
i-i'-Hubbard model combining the four-field slave-boson functional integral technique with the Bethe 
cluster method. The ground-state phase diagram evaluated self-consistently at the saddle-point and 
pair-approximation levels shows the experimentally observed suppression of magnetic long-range order 
in the favour of a paraphase with antiferromagnetic short-range order. In this phase the uniform static 
spin susceptibility consists of interrelated itinerant and local parts and increases upon doping up to the 
transition to the Pauli paraphase. Using realistic values of the Hubbard interaction we obtain the cusp 
position and the doping dependence of the zero-temperature susceptibility in reasonable agreement with 
experiments on Laa-^Sr^CuO,!. 



The spin correlations in high-T c superconductors 
probed by neutron and NMR experiments and by 
the spin susceptibility in the paraphase (show- 
ing, for La2_5Sr,5Cu04 (LSCO), a maximum in 
the doping and temperature dependences Q) are 
believed to be caused by a strong Coulomb in- 
teraction in the Cu02 planes which may result 
in a considerable antiferromagnetic short-range 
order (SRO). In our previous work || the role 
played by SRO in explaining the unconventional 
magnetic properties of the cuprates was investi- 
gated within the 2D one-band Hubbard model 
based on the scalar four-field slave-boson (SB) 
approach 

In this paper the SRO concept is futher elab- 
orated, where attention is paid to band-structure 
effects on the stability of magnetic SRO versus 
long-range order (LRO) and on the doping de- 
pendence of the spin susceptibility. In the SB rep- 
resentation the Hubbard model can be expressed 
as 1 
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The functional integral for the partition function 
is calculated in the radial gauge and the static 
approximation for the SB fields pi a , di, and the 

(2) 

Lagrange multiplier field X ia (enforcing the con- 
straint f\j ia = pt^ + djd.), where the e, fields 
are eliminated by the saddle-point approximation 
for Af ) (guaranteeing e\e t + d\d { + Eo-pLPwt = 



1) ||. Transforming away the SB fluctuations in 
the transfer term (proportional to the band nar- 
rowing z ia ), we get the free-energy functional 



(U d*d t - mvi + m^i) - Tr In 
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where nij = J2a a Piu an d n i ~ X " ' 
SB representations of the local magnetization and 
particle number, respectively, £, = J2a a ^Ha ^ s 
the internal magnetic field, and Vi — ^ CT X ia . 
In (^), h denotes the uniform external magnetic 
field. 

We incorporate the SRO beyond the PM sad- 
dle point by an expansion in terms of the local 
perturbation 
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G° ija = FT{[( W - v° + <r(f + h))/{z°Y - ej:]- 1 } 
is the PM Green propagator (the superscript "o" 
denotes the uniform paramagnetic (PM) saddle 
point), where 

= — 2i(cos k x + cos k y ) — At' cos k x cos k y (5) 

is the 2D tight-binding dispersion taking into ac- 
count the transfer integral between nearest (t) 
and next-nearest (t') neighbours. Now we treat 



the fluctuations of rrii = rhiSi, £j = £jSj (s, = ±) 
and of the charge degrees of freedom by the ansatz 
bi — * b Si with 6 £ {to, £, n, v,d= d*}. Transform- 
ing the functional (0) to an effective Ising model 
in the nearest-neighbour pair ((ij)) approxima- 
tion [Q we obtain 

(6) 



J SiSj , 



with 



a,cr— ± 



where the two-site fluctuation contribution 

^aa'a = < ^{ij)a ( a i : a j ) \ <*; = <* IS found to be 
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Im In 1 - G^Tj^a^G^Ti^ai) ,(8) 



and T ia = V ia (l — G° icr V^) -1 denotes the scat- 
tering matrix. Performing the Sj sum in the 
partition function with (||) we treat the SRO 
in the Bcthc cluster approximation. Thereafter, 
we determine the saddle point for all Bose fields 
b a G {to q , £ q , n a , u a , d a }. Correspondingly, in 
our theory the SRO is self-consistently described 
at the saddle-point and pair approximation lev- 
els at each interaction strength U and hole doping 
8 = 1 — n. 

In the h — limit (fh a = to) we obtain 
two possible paraphases ((sj) = 0): (i) the para- 
phase without SRO (PM; J = 0, m = 0) and (ii) 
the paraphase with antiferromagnetic SRO (SRO- 
PM; J < 0, to > 0). Let us stress that, at T = 0, 
the SRO-PM phase must be distinguished from 
the phase with antiferromagnetic LRO (denoted 
by AFM) having a finite sublattice magnetization 
m A = p A ^ — p A ^ = —m B which is determined 
from the A-B saddle point (|] and differs from 
the SRO amplitude to. 

In Fig. 1 the ground-state phase diagram 
is depicted, where only the phases describable 
by the scalar SB approach are shown and the 
tight-binding density of states is used (instead 
of a semielliptic one taken in Ref. 11). At large 
enough interaction strengths (U/t > 5 — 6) the an- 
tiferromagnetic LRO makes way to antiferromag- 
netic SRO in a wide doping region, where we ob- 
tain an AFM ^ SRO-PM phase transition of first 



order at S Cl and a SRO-PM ^ PM transition of 
second order at 8^. For U/t = 8 (being realistic 
for the cuprates (gj) and t' = we get 8 C1 = 0.04 
and S C2 — 0.26. The inclusion of the t' term in (^) 
(t'/t = —0.16 is realistic for LSCO) extends the 
stability region of the SRO-PM phase, since hop- 
ping processes along the lattice diagonals favour 
antiferromagnetic correlations . 




Figure 1: Ground-state phase diagram 

Within our theory, the suppression of mag- 
netic LRO at 6 C1 ~ 4% (U/t ~ 8) observed 
in LSCO may be related to the persistence of 
SRO in the paraphase. The SRO breaks down 
at S Cl ~ 30%. Note that just at this doping the 
superconductivity in LSCO disappears. Corre- 
spondingly, the unconventional behavior of high- 
T c cuprates may be due to the presence of mag- 
netic SRO which we suggest from our results to 
play also a role in the pairing mechanism. 

The uniform static spin susceptibility x nas 
to be calculated according to 



lim > 

h->o ^ 



dh 



dW a 
: dh 
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where m a = m a a, W a (h,h*, J) is the probabil- 
ity for the Ising spin projection a at the central 
site of the Bethe cluster, and h* is the effective 



Bcthc field. The first term in (||) describes the 
change of the magnetization amplitude with the 
applied magnetic field and gives mainly the 'itin- 
erant' contribution to \. The second term de- 
scribes directional fluctuations of the local mag- 
netizations ('local' contribution) and is finite only 
in the SRO-PM phase. Note that the 'itinerant' 
and 'local' properties are interrelated and deter- 
mine both contributions to the spin susceptibility. 

Figure 2 shows our results for the doping 
dependence of the zero-temperature susceptibil- 
ity X (S)- In the PM phase (S > 5 C2 ) the SB 
band-renormalized Pauli susceptibility has a pro- 
nounced doping dependence in two dimensions 
which is strongly affected by t'. In the SRO- 
PM phase (5 Cl < S < S C2 ), the Pauli suscepti- 
bility is suppressed due to the SRO-induced spin 
stiffness against the orientation of the local mag- 
netizations along the homogeneous external field. 
Accordingly, at S C2 a cusp in %(<J) appears, where 
the U / 1 and t' dependences of the peak-position 
emerge from Fig. 1. Since, for S Cl < S < S C2 , 
| J | decreases with increasing 6, the susceptibility 
increases upon doping. 
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Figure 2: Uniform static spin susceptibility at 
T = (taking the realistic value t — 0.3 eV). 



magnetic core (—9.9 x 10 -5 emu/mol) and Van 
Vleck (2.4 x 10 -5 emu/mol) contributions which 
can be taken as independent of doping and tem- 
perature over the limited parameter region stud- 
ied here Q|. The experimentally observed pro- 
nounced maximum at a hole doping of about 25% 
and the qualitative doping dependence of x are re- 
produced rather well by our theory without any 
fit procedure. 

From our results we conclude that the con- 
cept of magnetic SRO in strong-correlation mod- 
els may play the basic role in the explanation of 
many unconventional properties of high-T c com- 
pounds. As motivated by neutron scattering 
probing the correlation length over several lat- 
tice spacings, our theory should be extended by 
the description of a larger than nearest-neighbour 
ranged SRO. 

This work was performed under the aus- 
pices of Deutsche Forschungsgemeinschaft under 
project SF-HTSL-SRO. U.T. acknowledges the 
hospitality at the University Bayreuth. 
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In Fig. 2 we have also depicted the spin con- 
tribution to the magnetic susceptibility of LSCO 
at 50 K obtained from the experimental data on 
the total susceptibility by subtracting the dia- 



